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^ ; Abstract 

We considered a two dimensional three electron quantmn dot in a magnetic field in the Wigner 
. limit. A unitary coordinate transformation decouples the Hamiltonian (with Coulomb interac- 

tion between the electrons included) into a sum of three independent pair Hamiltonians. The 
eigen-solutions of the pair Hamiltonian provide a spectrum of pair states. Each pair state de- 
fines the distance of the two electrons involved in this state. In the ground state for given 
■ pair angular momentum m, this distance increases with increasing \m\. The pair states have 

• to be occupied under consideration of the Pauli exclusion principle, which differs from that for 
c/3 , one-electron states and depends on the total spin S and the total orbital angular momentum 

= ^ nii (sum over all pair angular momenta) . We have shown that the three electrons in 
the ground state of the Wigner molecule form an equilateral triangle (as might be expected) 

• only, if the state is a quartet {S = 3/2) and the orbital angular momentum is a magic quantum 
number {Ml = 3 m;m = integer). Otherwise the triangle in the ground state is isosceles. For 

. Ml = 3m -|- 1 one of the sides is longer and for Ml = 3m — 1 one of the sides is shorter than 

C ' the other two. 

o 

PACS numbers: 

73. 21. La , 73.63.Kv Quantum dots 
^ ■ 73.20.Qt Electron solids 

^ ! 31.30.Gs Jahn-Teller effect 
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I. INTRODUCTION 



Quantum dots are artificial atoms or molecules, where the electron num- 
ber, the scalar potential and the magnetic field are tunable and can provide 
favorable conditions for all kinds of fascinating effects (for recent reviews 
see E], [sl). One effect is the formation of Wigner molecules (WMs). 
In y it has been considered and quantitatively described in the simplest 
system, the two electron dot, using the exact analytical solutions for this 
system ji]. In the present paper, another effect, namely a Jahn- Teller-like 
distortion of the WM, is analyzed. Again analytical solutions for the sim- 
plest system, where this effect can occur (the three electron dot), proved 
useful. 

From the very beginning we should be aware of the fact, that Wigner 
localization in circular symmetric systems cannot be identified using the 
electron density n(r). Because the Hamiltonian commutes with the total 
angular momentum operator, the electron density can always be chosen cir- 
cular symmetric jB] and exhibits only a radial shell structure. Non-circular 
solutions are an indication of degeneracy, but not of Wigner localization. 
The fact that the electrons keep more or less fixed distances from each 
other due to a strong electron-electron correlation can be observed in the 
pair correlation function 
G(r) =<ij\j:6{r,-r,-r)\ij> 

i<j 

or the two-particle density matrix 
P(r,r') =< Y. S{r^ - r)6{rj - r')\ijj > 

i<j 

or in spin resolved versions thereof. For few electron dots both quantities 
are equivalent in exhibiting the existence of Wigner molecules. (In Q the 
pair correlation function G{r) in conjunction with the density n(r) was 
used for describing the effect instead of the nowadays favored two-particle 
density matrix P(r, r')). 

A useful and illustrative notion of a WM in a circular symmetric con- 
finement is a rotating and vibrating finite electron lattice |6|, 0,S,i. Due 
to this picture, a three-electron WM in a circular environment would form 
a equilateral triangle. In the present paper, however, we have shown that 
this is only the case in the classical limit where the distances between 
the electrons are so large that the exponentially decaying overlap between 
the localized wave functions of the individual electrons do not matter, but 
the long range Coulomb interaction is still effective. Shorter distances be- 
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tween the electrons and consequently overlap between the localized electron 
wave functions may lead to a distortion. Whether a distortion occurs de- 
pends on the angular momentum. This distortion of the seemingly natural 
equilateral symmetry is reminiscent of the Jahn-Teller effect, although in 
the present case the displaced objects are not atoms, but localized single 
electrons. (For a recent book on the Jahn-Teller effect including a lot of 



references see e.g. P^.) As shown below, the ground state of the WM 
for three electrons is isosceles including the equilateral as a special case. 
Excited states can be completely non-symmetric. 



Meanwhile there are a couple of papers p^, O, [12|, [13|, [15| which are 
focused on the case of three electrons. More work, which includes three 
electrons as a special case, can be found in recent reviews [El, The issue 
of the deformation of the WM without magnetic field has been discussed 
15| by investigating the two-particle density matrix P(r, r') using wave 



m 



functions from exact diagonalization in an oscillator eigenfunction basis. 
The author found that with increasing coupling parameter A = /o/as (/o = 
yjh/ {rrfujo), ttB = / {rrf e^)) there is a level crossing at A = 4.343, where 
the ground state switches from the {L,S) = (1, 1/2) to the (0,3/2) state. 
He developed the following qualitative picture. In the quartet state (ground 
state for large A), the WM forms a equidistant triangle for all A, whereas in 
the doublet state (ground state for small A) the spatial distribution is less 
trivial. Here, for infinite A the electrons occupy the edges of an equidistant 
triangle. With decreasing A the triangle is increasingly deformed into a 
isosceles one. This picture is a special case of the theory presented here. 
The purpose of the present paper is to extend the considerations to finite 
B and to present a simple quantitative model. 
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II. DECOUPLING INTO THREE PAIR PROBLEMS 



First, we want to review that part of the decoupling of the three-electron 
Hamiltonian into three pair Hamiltonians |ll3|, which is vital for the gen- 
eral understanding of the current paper. A numerical check of the validity 
of this decoupling procedure for the three-electron problem in a three- 
dimensional confinement without magnetic field can be found in |ll7| . Un- 



like in |ll3|, in the present paper we did not take advantage of the fact 
that for certain external field strength and quantum numbers there are 
analytical solutions of the pair equation, but instead we solved the radial 
pair equation (one-dimensional eigenvalue problem) numerically, whenever 
a concrete solution is required. The basic results, however, can be under- 
stood without having concrete numerical wave functions, but using only 
the Pauli principle. 

The Hamiltonian for three electrons in a homogeneous magnetic field 
B with the vector potential A(r) = (1/2) B x r and a harmonic scalar 
confinement (oscillator frequency uq) reads 

3 
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(1) 



i<k 

This is in atomic units h = m = = 1. For a model, where m is 
replaced by an effective electron mass m* and by the effective (screened) 
charge e*^ = e^/e, our results are in effective atomic units (a.w.*) defined 
hjh= rrf = e*^ = 1. Energies in our figures, which are generally given in 
units oi uq 01 Uc = B / c 01 a. combination thereof, are therefore independent 
of the background parameters. In order to avoid the dependence of the 
results on any material dependent parameters, the Zeeman term Hgpin = 

3 

{g* l^) /i^ cr, • B with /i^ = e*h/2Trfc, = ±1 is omitted. Besides, 

the Zeeman term has no influence on the focus of this paper, namely the 
spatial distribution of the electrons for given quantum numbers. It only 
shifts the energies and determines, what the quantum numbers of ground 
state are. 

Now, we apply a unitary transformation from the original position vec- 
tors Xi to new ones Xj 

1/3 a b 
b 1/3 a r2 (2) 
a b 1/3 







X2 
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where a = 1/3— 1/"\/3 and 6 = l/3+l/\/3. From the inverse transformation 
we obtain 

I"! - 1*2 = 

r2 - rg = \/3 ( X - xi ) (3) 




rs - ri = VS ^X - 



where X = ^X]i=i^i center of mass (cm.) in the new coordinates, 

which agrees with the cm. in the original coordinates R = ^ X^^^i ^i- We 
do not use the Jacobi transformation (see also Appendix C), which sepa- 
rates the cm. coordinate from the relative coordinates and which breaks 
the symmetry between the new quasi-particles. Instead, our transforma- 
tion retains the symmetry and the cm. is not an independent variable. 
Because any unitary transformation leaves the kinetic energy and the har- 
monic external potential invariant, the Hamiltonian in the new coordinates 
reads 

3 




(4) 



While being still exact, (Hj) is not completely decoupled because X contains 
all coordinates. In the Wigner limit, however, where the uncertainty of the 
cm. vector X is small compared with the mean electron- electron distance, 
we can neglect X in the denominator of the interaction term in ([11). As 
shown in the Appendices A and B, the quantitative errors introduced by 
this approximation are small and do not invalidate any of the qualitative 
conclusions of the present paper. Then, the Hamiltonian in zero order in 
X 

H^'^ = j2h^ (5) 
i=i 

decouples into a sum of three independent pair Hamiltonians, which can 
be rewritten as 

, = _lv^+1^2,2+l,^,^ + _L^ (6) 



Here, u = y^cjg + {uc/2y is an effective confinement frequency, Uc = B / c is 
the cyclotron frequency, and Iz is the orbital angular momentum operator. 
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This suggests the definition of a pair equation 

h (pq{x) = £q (^q(x) (7) 

with the normahzation condition J d^x |(/?g(x)p = 1. The subscript q com- 
prises aU quantum numbers. The pairs are essentially quasi-particles. 

In polar coordinates x = (x, a) we can make the following ansatz for the 
pair functions 



where m is an eigenvalue of Iz and therefore the angular momentum of a 
pair. Inserting (El) and (El) into ([71) provides the radial pair equation 



1 d 



2 



2 dx^ ' "^^^ 
with the effective pair potential 



u{x) = £ u{x) (9) 



^.//W = ip-y^ + i^V + ^ (10) 

and the definition 

1 , , 

£ = £ — -mujc. (11) 

The normalization condition reads J dx \u{x)\'^ = 1. Fig HI shows Veff{x) 

o 

for two typical effective confinement frequencies uj. Observe that Veff{x) 
for m = has a minimum at non-zero x only for uj < ujcr = "\/3/2 = 0.866. 

Because of the decoupling in zero order, the total eigenvalues and orbital 
eigenfunctions of H^^'> read 

^gi,g2,g3(xi,X2,X3) = (/^^^(xi) • (^52(x2) • ^q.i^s) (13) 

The total energy is a sum of pair energies and the total orbital eigenfunction 
is a product of pair functions. In the original coordinates, the orbital 
eigenfunctions read 

$,,,,,,3(ri,r2,r3) = 0,,(2 - 3) • 0,,(3 - 1) • ^,3(1 - 2) (14) 

where, for the sake of obtaining simpler formulae, we introduced the fol- 
lowing shorthand notation (j){i — k) = (p (li — (vi — rk)/VS). 
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FIG. 1: (Color online) Effective pair potential for two typical effective confinement frequencies 
Co. 



In order to avoid a basic misunderstanding, we want to stress that our 
independent-pair picture does not imply a division of all electrons into 
pairs, where each electron belongs exactly to one pair, as familiar from 



geminal approaches in Quantum Chemistry (see recent reviews in jl4j |) 
and which works only for even electron number. Instead, the total WF 
( fT4l ) is a product of two-electron (pair) functions, where each electron is 
involved in pairs with any of the other electrons. In |ll3| it has been shown 
that the high field limit of our wavefunction (WF) for R = agrees with 
the Laughlin WF for three electrons |jlo| • 

Now we want to discuss in an illustrative way, how the total WF ([HI) 
reflects the basic property of WMs, namely the strong e-e-correlation. Let's 
consider the limit R = X = first. The radial part of the pair functions 
u{x) shows a peak at a certain coordinate value around the minimum of 
Veff (see FigJSj), which we call the pair length and which will be calculated 
in Sect. 4. This means that the probability density |0gi(2 — 3)p of pair qi 
is peaked in r-space whenever |r2 — rs] equals the pair length. Considering 
the special cyclic structure of the total WF ( |T3l ) we can conclude: The 
probability density has peaks, if the three electrons form a triangle with 
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sides equaling the three pair length. The angular orientation of this triangle 
is arbitrary, and its translational location is defined by R = 0. This triangle 
is blurred by the finite width of the peak in u{x). 

Although our theory is aimed at the strong correlation limit, we want to 
point out that also the weak correlation limit can be described by product 
wave functions in the new coordinates. (This will also be demonstrated 
in the next Section by investigating the numerical solutions of the pair 
equation in this limit.) The reason for the decoupling of the Hamiltonian 
(IH), however, is quite different. Whereas in the strong correlation limit 
decoupling occurs because X can be neglected versus x^, this is not the 
case for weak correlation. In the weak correlation limit, however, the whole 
e-e-interaction term can be neglected versus the kinetic energy. This leads 
to decoupling as well. Describing non-interacting electrons by products of 
WF in the new coordinates rather than simply by one-electron functions 
is just a complicated (but equivalent) way of describing the same system. 
This ambiguity is a consequence of the fact that (after neglecting the e-e- 
interaction) our Hamiltonian is invariant under any unitary transformation 
of the coordinates. 

Because the pair functions are eigenfunctions of the orbital angular mo- 
mentum operator in x-space with eigenvalues m^, the total orbital eigen- 
functions (1131 ) are eigenfunctions of the total orbital angular momentum 
in x-space with eigenvalues Ml = Xli^i^*- Because the transformation 
back to the r-space is unitary, ( fMf ) has the same eigenvalues Ml like (fT3l). 
Moreover, through ([HI) and (fT^. Ml is linked to the parity of the total 
wave function. Even (odd) Ml means even (odd) parity, i.e. symmetry 
(antisymmetry) with respect to inversion of all coordinate vectors. 
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III. PAIR ENERGIES AND ELECTRON LOCALIZATION 



The pair energies £q in our pair approach are of the same central impor- 
tance for the electronic structure of the system as the one-particle energies 
in independent particle systems. The major qualitative difference between 
both approaches lies in the occupation (Pauli principle) of the energy lev- 
els, which has been investigated in Ref.|13[. Figs J21I11 show the dependence 
of the ground state energies of the pair levels on the orbital angular mo- 
mentum m of the pairs for fixed external field. 

FigO applies to the limiting case of vanishing magnetic field. In the 
limit cjQ ^ oo the pair levels converge to the one-electron levels for non- 
interacting electrons (Fock-Darwin levels) 

^non-^nt = (2 n + |m| + 1) + LJc (15) 

where the radial quantum number n = 0,l,2,-- - is the degree of excitation 
for given angular momentum m. This can be seen in the radial pair equa- 
tion (I91JT01): For large uj the state is compressed into the region of small x, 
where the centrifugal term (oc l/x^) is much larger than the e-e- interaction 
term (oc 1/x) and the latter can be neglected. Without the e-e-interaction 
term the pair equation agrees with the one-electron Schrodinger equation. 

FigJ21 shows the complementary case of vanishing confinement. In the 
limit of strong magnetic fields, these pair levels approach the ground states 
(for given m) of the Fock-Darwin levels (ITHI). The analytical verification of 
this statement is evident in eqs. (l9lJT0l). It is also clear intuitively, because 
the strong magnetic field out-plays the effect of the Coulomb interaction. 
A case with a fixed finite confinement cjq, which is of the order of magnitude 
of quantum dots in GaAs (in effective atomic units) is shown in FigJH The 
angular momentum rrimin with minimum pair energy is now non-zero. 

Fig J5] provides the angular momemtum rrimin of the lowest pair energy for 
given magnetic field and parabolic scalar confinement. The lines separate 
the regions with adjacent rrimin- In the upper left part of the plot mmin = 
and and in the lower right part rrimin converges successively to — oo. In 
between only the lines separating the 10 lowest rrimin are shown. Although 
only a limited region of the parameter space is shown, the general features 
are obvious: All phase boundaries in the log-log-plot are well represented 
by parallel straight lines log{iOQ) = log{A) + B ■ log{LOc). Consequently, the 
boundaries in a linear plot are power functions uq = A-u^ where B ~ 0.75 
is universal and the factor A depends on the boundary in question. 
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FIG. 2: (Color online) Pair energies (lowest state for given m) for zero magnetic field (wc = 0) 
for the three confinement frequencies indicated on the levels. The colored values are obtained 
by numerical solution of ([9j). The thin (black) levels show the one-particle energies for non- 
interacting electrons given in (|15p . The lines connecting the levels are a guide for the eye. 

FigJSl shows how electron locahzation can be visuahzed in our approach. 
Locahzation means that (the radial part of) the pair function which 
describes the electron- electron distance, is peaked at a finite x. There are 
two reasons for localization: 

a) for small coq (and small or vanishing Uc) the electrons are pushed away 
from each other by the (last) interaction term in ([TUI). On the other hand, 

b) for large cUc (and small loq) the modulus of the angular momentum of the 
ground state is large (see FigJS)- In this case the separation of the electrons 
is caused by the (first) centrifugal term in (fTUf). This mechanism does not 
work for vanishing ujq because then the angular momentum is ill defined. 
It does not work for vanishing electron- electron interaction either. 
Therefore, electron-electron interaction alone can localize the electrons 
without the assistance of the magnetic field, but the magnetic field alone 
cannot do the job. As seen in FigJSl the strongest localization is gained for 
small a;o(=0.01) and large (or at least medium) Uc. On the other hand, a 
medium cjc{=l) cannot achieve anything if cjq is of comparable size. (The 
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FIG. 3: (Color online) Pair energies (lowest state for given m) without confinement potential 
{ujQ = 0) for the three cyclotron frequencies indicated on the levels. The thin (black) levels show 
the one-particle energies for non-interacting electrons. 



curve for a;o=l, c<;c=0.01 agrees almost completely with the curve for a;o=l, 
ujc=l.) In this case we need a very strong magnetic field a;c(=100) for lo- 
calizing the electrons. A quantitative measure for localization is the mean 
square deviation of x from the expectation value < x >. 

/•oo 

{Ax)^ = dx [u{x)f {x- < X >)^ 
Jo 

The results read in the same order as shown in FigJS {Ax^ = 0.0749, 
0.00155, 0.1987, 0.2017, 0.00717. 

IV. PAIR LENGTH AND DISTORTION OF THE WIGNER MOLECULE 

As seen in FigJTl the pair length define the distortion of the Wigner 
molecule. The pair length, or mean electron-electron distance, can be 
obtained from the following quantities: 

(i) the minimum position of the effective pair potential Veff{x) 

(ii) the maximum position of the radial pair function u{x) 
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FIG. 4: (Color online) Pair energies (lowest state for given m) for a dot confinement (wq = 0.2) 
for the four cyclotron frequencies indicated on the levels. 

(iii) the average electron distance 



< X >= 




The real electron- electron distance Ar is obtained from the above defined 
values by multiplication with as seen from ([31). Definition (iii) agreed 
with (ii) if the radial pair function was symmetric with respect to the 
maximum, and (ii) agreed with (i) if it was (5-function like. In our curves 
discussed below, the radial wave functions of the ground state for the 
corresponding m have been used throughout. Observe that for m = and 
uj > -\/3/2 the first definition breaks down completely, because Veff{x) has 
no minimum at non-zero x (see also FigHj). Fig JHl compares these definitions 
for 171 = 0. As to the agreement of the first and simplest definition with the 
two more sophisticated ones one has to consider, that a discrepancy of this 
size occurs only for m = 0, where the centrifugal potential in 2D is negative. 
The difference is much less for \m\ > 0. Fig J21 shows the dependence of the 
pair length Ar on the orbital angular momentum m. Apart from getting 
an idea about the order of magnitude of the dependence, we notice that 
the pair length is the larger the larger \m\ is. (Observe that Veff{x) and 
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FIG. 5: Orbital angular momentum rrimin of the pair with minimum pair energy for given 
cyclotron frequency oJc and confinement frequency ujq. 

u{x) depend only on the modulus of m.) This follows directly from the 
fact, that the centrifugal potential in (fTUj). which grows with growing |m|, 
pushes the electrons within a pair away from each other. This effect can 
also be understood on the basis of classical mechanics. From the almost 
linear behavior in the log-log-scale we conclude, that the function of Ar 
versus uj is roughly a power law. 

For ground states, the qualitative features of the distortion may be 
grouped into three cases. In accordance with Fig.s EHU we assume for 
the lowest pair state rrimin < 0. 

(a) If all three pairs are identical, i.e. Ml = 3 rrimin: the triagle is equilat- 
eral. This case can occur only in quartet states (parallel spins) and these 
total orbital angular momenta are called 'magic'. 

(b) The thick (black) lines in FigJTUl describe the angular distortion of a 
WM composed of two pairs with iJimin = and one pair with m = — 1. 
Qualitatively, all pictures for two pairs with angular momentum mmin and 
one with rrimin — 1, i-e. with Ml = 3 rrimin — 1, agree. Then the WM spans 
a triangle where two sides are equal and one side is longer than the other 
two, because \mmin ~ M > |^mm| for rrimin ^ 0. The Pauli principle for 
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FIG. 6: Localization of the pair function for the external fields given in the caption, mmin is 
the angular momentum of the state with lowest energy and < x > is the expectation value 
< X >= J dx X u{x)'^. The curves are renormalized to conserve the normalization. 



pair states |ll3| demands that this case can occur as in quartet (non-magic 
angular momentum) as well as in doublet states. 

(c) FigHni also shows the case for two pairs with rrimin = — 1 and one pair 
with m = 0. All WMs with two pairs with rrimin < and one pair with 
^?7im + 1 providing Ml = 3 rrimin + 1 look similar. One side of the trian- 
gle is shorter than the other two equal sides because Im^m + 1| < Im^i^l 
and because this case can happen only for rrimin < 0. The Pauli principle 
imposes the same restrictions to the quantum numbers as in case (b). 

It is clear that the distortion has to vanish in the limit cj ^ 0, because 
in this limit the electrons behave like classical particles and their localized 
wave functions do not overlap. The distortion, however, is a quantum 
mechanical effect. 

As can bee seen from the Pauli principle and the minimization of the 
total energy, in ground states the three pairs are never all different from 
each other. This does not apply to excited states where all three angular 
momenta can be different leading to a completely non-symmetric WM. The 
results of our theory for excited states have to be considered with some 
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FIG. 7: (Color online) Schematic picture of the distorted three-electron Wigner molecule with 
total orbital angular momentum Ml = —1 (left) and Ml = —2 (right). The thin lines depict 
the undistorted WM. 




I > > > > > > > > > 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

0) 



FIG. 8: (Color online) Pair length for m = as estimated from different procedures (see text). 

caution, however, because its applicability to excited states (in particular 
cm. excitations) has not yet been investigated thoroughly. But this is not 
the topic of this paper. 
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log CO 

FIG. 9: (Color online) Pair length calculated from < x > for several pair orbital angular 
momenta \m\. 

V. SUMMARY 

Using a simple unitary coordinate transformation and the Pauli princi- 
ple, we have shown that a three-electron Wigner molecule in a quantum dot 
and a magnetic field can show a Jahn-Teller-like distortion. The qualita- 
tive fact of distortion can be shown analytically without any computations, 
but for a quantitative estimate we have to solve a one-dimensional eigen- 
value problem (radial pair equation) numerically. In the ground states the 
Wigner molecule is either equilateral or isosceles, whereas excited states 
include completely non-symmetric geometries as well. If the state is a 
quartet = 3/2) and the orbital angular momentum is a magic quantum 
number [Ml = 3 m;m = integer), the triangle is equilateral. In the dou- 
blet state ('S' = 1/2) and for = 3m + 1 one of the sides is longer, and 
for Ml = 3m — 1 one of the sides is shorter than the other two. 



16 




FIG. 10: (color online) Angles in the distorted Wigner molecule calculated from < x >. The 
thick (black) lines refer to the case where two pairs have m = and one pair has m = —1 as 
sketched in FigjT] (left). The thin (red) lines belong to the case where two pairs have m = — 1 
and one pair has m = as shown in FigjT] (right). As in FiglTJ a gives the values of the two 
equal angles and (3 is the third angle. 



VI. APPENDIX 

A. Estimate of the ratio x/R 

For B = the ratio x/R, which is the relevant parameter for the decou- 
phng in ([Ij), can be estimated easily. The definition of the new coordinate 
X in dSj) reads 

x = R-^Ar (17) 

where we dropped the indexes and denote the electron distance by Ar. An 
idea of the e-e-distance can be obtained from the classical e-e-distance in 
the ground state which is Avd = 3^/^ Uq The width Rq of the cm. 
probability distribution can be deduced from the exactly known cm. wave 
function (see below) providing 



Rn = J d^n R |$c.m.(R)|' = (v^/6) Uq'^^ 



Formula ( fTTj ) with Ar/"\/3 > R provides for the modulus 

1 Avci X < ^ < ^ ^^^^ + 1 (18) 

Rq ~ R ~ Rq 
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Because Ar,;/i?o = (6 • Z^l'^jix^l^) cjq ^^'^ ^ (X) in the Wigner limit cjq ^ 0, 
we conclude x/ oo and for the small decoupling parameter R/x 0. 

B. Validity of the approximation X=R=0 

Here we want to estimate the error due to the neglect of X = R in 

i) the e-e-interaction term in the transformed Hamiltonian (jll) 

ii) the pair length considered in Sect. 4. 



In a previous paper the e-e-interaction term was expanded in a multi- 
pole series (for X < x^) and the correction terms were considered in per- 
turbation theory. The problem with this approach is that the resulting 
series converges slowly. Additionally, it provides only the total energy in a 
simple way. Therefore, we adopted here a different approach. Because the 
probability distribution w{X) of the cm. X = R is known exactly from 
Kohn's theorem (see Appendix C) we averaged the corresponding quanti- 
ties containing X = R with w{X) as weight function. (For the estimates in 
this Appendix we adopted for the angular momentum of the cm. system 
m = 0.) 

We replaced the e-e-interaction potential 
by the averaged potential 

Fee(x,) = j d^X W{X) ■ yee(x„ X) (20) 

This virtually means using a smoothed e-e-potential (see FigHT]). Unlike 
perturbation theory, this approach does not destroy the independent-pair 
picture. Putting R = X = can also be viewed as using a (5-function-like 
cm. distribution as found in the ground state of classical mechanics. 

Despite the fact that the averaged e-e-potential in FigHT] deviates or 
small X strongly from the curve for R = 0, the change in the radial part 
of the pair function is very small. This comes from the fact that the 
probability density [w(x)]^ is small for those x, where the averaged e-e- 
potential is changed by averaging. For |m| = 5 the two curves for R = 
and w{R) are hardly distinguishable. The maximum shift in the pair 
energies due to averaging shown in Fig|I2]is only about 1%. 

Now we are considering the pair length. The relation between electron 
distances — r^, the new coordinates x/ and the cm. vector R is given in 
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FIG. 11: (Color online) Comparison of the e-e-interaction part of the effective pair potential for 
R = (thick, full, black) with the result using the cm. distribution w{R) (thick, broken, black) 
for ui = 0.1. The radial parts of the pair function ([9]) for \m\ = and 5 and for R = (full) and 
w{R) (broken) are colored. 



Denoting — for a chosen pair by r omitting the index at x;, we 
have 

r' = -i=r = R-x (21) 

Under the assumption of statistical independence, the probabihty density 
Wr{r') foUows from the known probabihty densities 

and 

wn{R) = ^[^(i^)]' (23) 

using 

Wr{T') = j d^R ujr{R) w;,(|R - r'l) (24) 

For the definition of the pair length we use the expectation value of the 
e-e-distance 

<r' >= [ Sy' r' Wr{r') (25) 
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FIG. 12: (Color online) Comparison of the pair energies (lowest state for given m) for for R = 
(colored) and for the cm. distribution w{R) (black). The confinement frequency is the same as 
in FigUJand the cyclotron frequencies are indicated at the levels. 



Using (EID, rescaling < r >= < r' > and after changing one of the 
integration variables we obtain 

< r >= a/S J (fyi x{x) Wx{x) = Vs j dx x{x) [u{x)]^ (26) 

with the weight function 

j d^R wr{R) |x - R| (27) 



For 5-function-like cm. distribution (corresponding to the approximation 
R = 0) we have x{x) = x and the result from (ESI) agrees with the result 
from (fT6l). Fig]T3l shows that x{x) calculated with a finite distribution for 
R deviates only for small x from the result with R = 0. In this region, 
however, the radial pair function is small and the pair length for both 
approaches shown in Fig|8] deviate only marginally. 

In the end we want to emphasize that the main conclusions of this pa- 
per, namely the qualitative statements about the distortion of the WM, 
are not influenced by averaging over R. There are only small shifts in 
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FIG. 13: (Color online) Weight function x{x) for the calculation of the pair length for R = and 
for the finite distribution w{R) (black) as well as the radial pair functions with R = (colored) 
for the pair momenta |m| = and 5. 



the pair energies and pair length. There is a hand-waving argument for 
this robustness against changes in the cm. vector. It is connected to the 
Generalized Kohn theorem, which states that the cm. coordinate decou- 
ples exactly from properly defined relative coordinates (see Appendix C). 
Therefore, no matter how large the width of the cm. distribution is, it has 
no influence on the corresponding relative coordinates and on the internal 
structure of the WM. However, we have to consider that our relative co- 
ordinates are not completely decoupled from the cm. vector, but are 
coupled weakly as indicated by the small shifts. 



Jacobi transformation 



Applying the Jacobi transformation (see e.g. |19|) 

1 -1 
1/2 1/2 -1 
1/3 1/3 1/3 



Pi' 








R 





(28) 
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to ([ID decouples the cm. vector R from the relative coordinates 



H = H,m(R, P) + Hrelip,, TT,) (29) 

where P = ^ - Pi is the total momentum and tTj are the canonical momenta 
belonging to p^. The distance between the electrons and thus the internal 
structure of the WM depends only on the relative coordinates p^ 

ri - r2 = Pi 

1 , , 

r2 -ra = Pi + P2 (30) 

1 

rg-ri = --P1-P2 

This means that there is no correlation between the cm. vector and the 
structure of the WM. For the estimates in Appendix A we need only the 
cm. Hamiltonian which reads (for arbitrary electron number N) 



H -1 



^P2 + 1 cD^ (iVR)2 + i cu, (iVR) X P 



(31) 



The cm. eigenvalues are independent of N and agree with (fT5|). but the 
eigenfunctions are homogeneously compressed by a factor of N . The latter 
read in polar coordinates R = (i?, a) 



$„,(R) = -=7^H(i?) (32) 

V Z7T 



with the radial part for the lowest state for given m 



y |m|! 

With ( 1321 ) and ( I33l ) the probability density for R reads 

Wm{R) = (l/27r) [n\^\{R)f. (34) 
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